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ABSTRACT 


An  analytical  approach  to  particle  erosion  of  solid 
surfaces  Is  developed  based  on  a  model  of  erosion  pit 
nucleatlon  and  growth.  The  Hertzian  theory  of  impact  Is 
used  to  Introduce  the  kinematic  and  material  parameters 
Into  this  general  statistical  formulation  of  erosion 
processes.  In  order  to  minimize  the  number  of  unknown 
factors  In  the  Initial  development  of  a  predictive  erosion 
theory,  the  particular  case  of  solid  spherical  particles 
Impacting  directly  on  the  plane  surface  of  an  elastic 
half-space  Is  Investigated. 

Erosion  tests  were  conducted  on  the  AFML-Bell  erosion 
apparatus  in  conjunction  with  the  modeling  studies.  The 
erosion  mechanisms  In  fused  silica  have  been  Identified 
when  eroded  by  multiple  Impacts  of  spherical  glr^ss  beads. 
The  erosion  process  for  2.8  mil  beads  Impacting  at  200  fps 
was  found  to  differ  from  the  sequence  of  erosion  mechanisms 
associated  with  11  mil  beads  at  the  same  velocity.  The 
observed  erosion  mechanisms  could  be  correlated  with  the 
measured  erosion  rates. 

The  erosion  of  fused  silica  can  be  modeled  In  terms 
of  a  preferential  pit  nucleatlon  and  growth  process.  The 
quantitative  evaluation  of  the  statistical  parameters  In 
the  erosion  model  Is  being  Investigated.  The  controlled 
erosion  tests  provide  evidence  for  the  applicability  of 
the  Hertzian  relations  to  the  Impact  conditions  of  Interest 
In  solid  particle  erosion. 
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I.  INTRODUCTION 


The  erosion  of  engineering  materials  moving  through  rain, 
sand,  and  dust  environments  poses  serious  operational  and  main¬ 
tenance  problems  In  a  variety  of  military  and  commerlcal  systems. 

To  counter  this  problem,  the  Air  Force  Materials  Laboratory  has 
established  programs  of  erosion  research  to  develop  advanced 
erosion  test  apparatus,  to  characterize  the  erosion  of  aerospace 
materials  In  known  environments,  to  establish  the  basic  mechanics 
of  erosion  engineering  materials,  and  to  develop  new  aerospace 
materials  having  Improved  erosion  resistance. 

The  complexity  of  the  loading  conditions  associated  with 
erosive  environments  and  the  lack  of  knowledge  of  the  general 
behavior  of  material  surfaces  subjected  to  these  loading  conditions 
have  been  major  obstacles  to  the  detailed  analysis  of  erosion 
phenomena.  In  order  to  confront  this  pacing  problem,  attention 
In  this  report  has  been  focused  on  the  simplest  combination  of 
circumstances  which  lead  to  a  realistic  model  of  the  erosion 
process  for  particle  Impacts. 

The  development  of  mechanistically  sound  erosion  models 
complements  the  overall  Air  Force  erosion  research  effort  In  two 
ways.  First,  such  models,  even  If  only  approximate,  will  help 


significantly  In  the  extrapolation  of  erosion  test  results  to 
service  conditions,  and  hence  facilitate  the  engineering  choice 
of  materials  for  use  In  particular  service  environments.  Second, 
the  development  of  mechanistically  sound  erosion  models  will 
allow  estimations  to  be  made  of  the  relative  Importance  of  the 
various  particle,  material,  and  Interaction  parameters  which 
Influence  the  complex,  cumulative  process  of  erosion,  and  hence 
aid  In  the  Interpretation  of  erosion  mechanism  studies  and  In 
the  design  of  new  materials  or  coatings  having  Improved  erosion 
resistance  In  specific  service  environments. 
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II.  SUMMARY 

.  ,  I  I  '  . 

■ 

I  ^ 

‘  •  I  , 

The  erosion  of  brittle  materials  by  multiple',  solid  ; 
particle  Impacts  has  been  Idealized  as  a  process  of  pit  nucleatlon 
and  growth  on  the  surface  of  the  eroding  solid.  An  analytical 

I 

model  Is  developed  based  on  the  assumption  that  the  erosion  pits ' 

I 

J 

are  rather  uniform  In  depth  during  the  initial  stages  of  growth 

i 

and  that  typical  pits  undergo  appreciable  lateral  growth  before 
conjoining  with  other  growing  pits.  'The  eroding' solid  Is  cop- 

I 

sldered  to  be  a  laminate  of  plates  whose  Individual  thickness  Is 
determined  by  the  depth  of  propagation  of  i  typical  Impact  fracl^ures. 

I 

The  governing  equations  have  been  programmed  for  digital  computer' 
computations  of  the  erosion  'rates  and  weight  loss  as  functions 
of  exposure  time. 

>  I  '  . 

The  mathematical  formulation  requires  the  specification  of  ' 
a  time-dependent  pit  nucleatlon  probability  fiiinctlon  and  growth 

I 

rate  for  each  layer  of  the  laminated  model  of  the  eroding  sdlid., 

The  value  of  the  theory  as  a  predictive  tool,  In  the  development 
of  erosion  resistant  materials  depends  on  how  much  Information 
pertaining  to  the  kinematics  of  the  eroding  environment  hnd  the  . 
materials  characteristics  of  the  eroding  solid  can  he  Incorporated 
into  the  explicit  analytic,  representation  of  the  pit' nucleatlon 
and  growth  rate  functions.  Ihe  direct  Impact  of  small,  spherical 

•I  I  ‘ 


I 


} 


glass  beads  on  the  plane  face  of  relatively  flaw-free  fused 
silica  (Corning  79^0)  samples  was  Investigated  In  order  to 
minimize  the  number  of  unknown  factors  pertaining  to  the  kine¬ 
matic  and  material  parameters.  The  Impact  conditions  were 
selected  so  that  results  from  the  Hertzian  theory  for  the  Impact 
of  a  spherical  elastic  body  on  an  elastic  half-space  could  be 
used. 

The  pit  nucleatlon  concept  Is  dependent  on  the  manner  In 
which  the  material's  surface  Is  fractured  by  the  multiple  solid 
'  particle  Impacts.  The  nucleatlon  of  an  errs  Ion  pit  may  occur 
through  one  of  several  processes  whose  relative  probabilities 
are  difficult  to  access  without  evaluating  a  detailed  model. 

This  Information  was  obtained  through  microscopic  observations  of 
the  erosion  mechanisms  In  a  series  of  controlled  experiments  using 
the  AFML-Bell  erosion  apparatus. 

Three  erosion  tests  were  conducted  on  fused  silica  saunples: 

2.8  mil  beads  Impacting  at  200  fps;  11  mil  beads  Impacting  at  200  fps; 
and  11  mil  beads  Impacting  at  300  fps.  The  experimental  program 
was  designed  to  obtain  Information  on  the  erosion  mechanisms, 
to  obtain  explicit  forms  of  the  pit  nucleatlon  and  growth  rate 
,  functions,  and  to  obtain  weight  loss  data  for  c'mp«u*lson  with  the 
predictions  of  the  theory. 

i  I 

I 


I 
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Entirely  different  eroalon  mechanisms  were  found  for  the 
2.8  and  11  mil  beads.  The  nucleatlon  of  pits  on  the  surface  of 
the  fused  silica  specimen  Impacted  by  2.8  mil  beads  was  due  to 
the  production  and  Intersection  of  ring  fractures  which  were  on 
the  order  of  one  mil  In  diameter.  Material  was  removed  from  the 
surface  layer  to  a  maximum  depth  of  10  microns  by  the  formation  of 
Irregularly- shaped  pits.  The  growth  of  the  nucleated  areas  was 
preferential  In  accordance  with  the  erosion  model. 

When  the  11  mil  beads  were  used,  there  was  a  much  higher 
probability  for  the  production  of  conical  fracture  subsurfaces. 

The  average  diameter  of  the  Inner  ring  fracture  seen  on  the  surface 
of  the  specimen  was  about  2.5  mils.  As  the  number  of  damaging 
Impacts  Increased,  the  subsurface  fractures  became  more  obvious. 

No  new  ring  fractures  were  formed  In  the  Interstitial  regions 
between  the  ring  fractures  produced  during  the  early  stages  of 
the  erosion  process.  Instead  a  very  dense  network  of  Internal 
fractures  developed  with  very  little  material  removal.  At  a 
certain  point  In  the  development  of  these  fractured  regions,  mater¬ 
ial  was  removed  to  a  depth  of  60  microns  for  Impacts  at  200  fps 
and  90  microns  for  Impacts  at  300  fps.  This  was  a  nucleatlon  and 
growth  process,  although  the  regions  observed  were  quite  erratic 
In  their  growth  patterns.  The  Interstitial  material  was  removed 
at  a  higher  erosion  rate  than  the  one  prevailing  before  this 
began.  After  most  of  the  Interstitial  material  was  removed,  the 
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eroded  surface  consisted  of  a  continuous  array  of  conical  frustums 
protruding  from  the  bulk  material.  The  majority  of  the  plateau 
regions  were  still  unaffected  by  the  particle  Impacts  and  the 
side  surfaces  were  quite  smooth.  At  this  point,  material  was 
removed  at  a  uniform  rate  by  the  fine- scale  chipping  of  the  plateau 
regions  of  the  frustum  population.  The  time  of  occurrence  of  the 
observed  erosion  mechanisms  can  be  correlated  with  the  weight 
loss  data. 

The  diameters  of  the  ring  fractures  In  all  three  tests 
were  found  to  be  In  general  agreement  with  the  Hertzian  predictions 
for  the  Impact  conditions  used  In  these  experiments.  The  develop¬ 
ment  of  a  suitable  fracture  criterion  will  require  further  study; 
however,  on  the  basis  of  the  results  obtained  here,  the  energy 
balance  approach  Is  preferable  to  the  statistical  flaw  theories. 

An  experimental  procedure  was  devised  for  direct  evaluation 
of  the  pit  nucleatlon  and  growth  rate  functions.  In  the  case 
of  the  2.8  mil  beads,  a  number  of  ring  fractures  with  a  diameter 
of  3  mils  were  observed  on  the  surface  of  the  sample  and  were 
Judged  to  be  the  major  contributor  to  material  removal.  The  fine- 
scale  chipping  due  to  the  smaller  ring  fractures  was  considered 
to  be  of  secondary  importance.  At  an  Intermediate  stage  of  the 
erosion  process  the  frequency  of  occurrence  of  the  larger  rings 
approached  zero.  It  was  then  discovered  that  the  fine- scale 
chipping  was  the  principal  erosion  mechanism.  There  was  no  growth 


-6- 


of  the  eroded  regions  associated  with  the  3  mil  rings.  No  useful 
data  were  obtained  pertaining  to  the  evaluation  of  the  statistical 
parameters  for  this  test.  The  reason  for  the  occurrence  of  the 
anomalous  3  mil  ring  fractures  la  being  Investigated. 

The  erosion  behavior  of  the  multiple  Impacts  of  the  11  mil 
beads  was  somewhat  unexpected.  The  experimental  procedure  for 
evaluating  the  statistical  parameters  was  not  suitable  for  this 
purpose.  Alternative  procedures  for  obtaining  the  necessary  data 
for  the  model  calculations  are  being  considered. 
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III.  DEVELOPMENT  OF  GENERAL  EROSION  THEORY 


An  analytical  model  for  describing  erosion  processes  In 
materials  which  erode  by  the  nucleatlon  and  growth  of  pits  In 
the  material  surface  is  currently  under  development.  The  eroding 
solid  Is  represented  by  a  laminate  of  Individual  plates  as  showi^ 
In  Fig.  1.  The  thickness  of  each  plate  Is  determined  by  the  meam 
pit  depth  for  the  pits  which  are  formed  at  each  stratum  In  the 
solid.  Assume  erosion  Initiates  through  the  nucleatlon  and 
growth  of  flat-bottomed  pits  which  penetrate  the  topmost  plate 
(Fig,  1).  Erosion  of  the  next  lower  plate  proceeds  through  the 
nucleatlon  and  growth  of  similar  pits  over  that  portion  of  the 
plate  surface  which  has  been  exposed  to  erosive  Impact  through 
regression  of  the  upper  plate  (Fig.  1).  Erosion  of  this  plate 
exposes  a  portion  of  the  surface  of  the  next  lower  plate,  which 
erodes  In  the  same  way  to  expose  a  portion  of  the  surface  of  the 
next  lower  plate,  and  so  on. 

Consider  first  the  erosion  of  the  topmost  lamella  of  the 
material.  Let  the  thickness  of  the  plate  be  D  and  Its  surface 
area  be  A.  Let  the  surface  be  exposed  to  erosion  at  time  zero. 
Presume  that  the  probability  that  an  erosion  pit  will  be  formed 
In  an  arbitrary  unit  area  of  surface  In  the  time  Interval  between 
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(a)  Solid  imaginad  at  a  Laminata  of  Platts 


(b)  Idtaliatd  Pit  Nucitation 


( 


(c)  Idtalitad  Pit  Growth 


(d|  Htaaiiiad  Mtchanitm  of  Pit  Dtaptning 


Fifure  I .  An  Idealization  of  Erosion  through  the  Nucication  and  Growth 
of  Pits  in  Successive  Layers  of  an  Eroding  Solid 


the  times  t  and  t  +  dx  after  Initial  exposure  is  known  and  is 
given  by  the  expression  Ii(T)dT.  The  function  Ii(t)  is  restricted 
in  that  it  must  be  integrable  and  must  vanish  in  the  domain  i  0. 
Further  assume  that  the  pit  nuclei  have  negligible  volumes  and 
they  grow  uniformly  over  the  surface  at  a  constant  rate  0..  remov¬ 
ing  material  to  depth  D  with  each  Increment  of  growth. 

If  erosion  pits  are  small  and  so  dispersed  over  the  sur¬ 
face  that  they  do  not  impinge  on  one  another,  then  a  pit  nucleated 
at  time  x  will  have  grown  to  cause  the  removal  of  a  volume  6v  at 
time  t. 

6vi(t,x)  ■  Daf(t-x)*  (1) 

where  Gi  is  the  lateral  growth  rate  for  the  first  layer.  It  is 
assumed  that  the  eroded  area  at  the  bottom  of  a  nucleated  erosion 
pit  at  depth  D  is  the  same  as  at  the  surface.  It  is  unlikely  that 
this  will  be  the  case  In  actual  erosion  processes.  The  Influence 
of  the  shape  of  a  nucleated  erosion  pit  on  the  volume  of  material 
removed  will  be  taken  into  account  after  microscopic  observations 
of  eroded  surfaces  indicate  the  form  of  the  shape  factor  adjust¬ 
ment  required.  The  volume  removed  at  time  t  due  to  the  growth  of 
all  pits  nucleated  in  the  time  Interval  (x,x+dx)  is  then 

6V|(t,x)  ■  AIi(x)6vi (t,x)dx  (2) 

However,  after  significant  erosion  takes  place,  the  uneroded  sur¬ 
face  area  is  much  less  than  the  total  surface  area  A,  and  so  the 
number  of  new  pits  nucleated  is  much  smaller  than  the  number  given 
by  AI|(x)dx.  Also  growing  pits  may  coalesce  with  one  another  to 
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form  single  larger  pits.  To  account  for  these  effects,  we  use 
the  net  uneroded  area  at  any  Instant  t  In  calculating  the  number 
of  pits  nucleated  ^n  the  time  Interval  (TiX+dt).  Hence,  Eq.  (2) 
becomes , 

6V,(t,T)  -  Gl(t-T)»Ii(T)(V.AV,(T))  (3) 

where  AV|(t)  ■  DAAi(t)  Is  the  total  voluire  eroded  from  the  start 
of  the  erosion  process  to  time  x  and  AAi(t)  Is  the  eroded  surface 
area  associated  with  the  same  time  Interval.  At  present  AVi(t) 

Is  an  unknown  quantity.  V  ■  AD  Is  the  volume  of  the  first  plate. 
Integrating  Eq.  (3)  from  (0,t), 

ft 

AVi(t)  •  Q\  (t-T)*Il(T)[V-AVi(T)]dT  ,  (k) 

0 

yields  a  Volterra  Integral  equation  of  the  second  kind  for  deter¬ 
mining  the  total  eroded  volume  of  the  first  lo'or  at  time  t  due 
to  all  the  pits  nucleated  during  the  total  time  of  exposure  to 
the  eroding  environment. 

We  follow  a  similar  procedure  for  evaluating  the  volume 
eroded  from  the  second  layer,  except  that  the  surface  of  layers 
below  the  first  Is  only  gradually  exposed  to  the  erosive  environ¬ 
ment.  For  convenience  of  presentation,  let  us  assume  that  the 
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mean  pit  depth  Is  also  D  for  the  second  layer.  In  accordance 
with  earlier  considerations  the  net  exposed  area  on  the  surface 
of  the  second  layer  at  any  time  ^  Is  AAi(#)  -  where 

AVi(<i)  ■  DAAi(#)  Is  the  unknown  quantity  which  must  be  determined. 
The  amount  of  new  area  which  Is  exposed  during  a  time  Interval 
can  be  expressed  In  terms  of  a  Taylor's  expansion, 

[AAi(<»4d*)  -  AAa(*+d*)]  -  [AA  i  (  ♦>-AAt  (♦)  ] 

dAAt(t(>)  dAAa(<r) 

:  — —  d*  -  — —  dt  (5) 

The  number  of  pits  nucleated  at  time  t  due  to  new  area  exposed 
during  the  time  Interval  Is  given  by 

/'dAAi(*)  dAAa(t)\ 

6n(T,t)  ■  - J  (6) 


It  should  be  noted  that  a  different  pit  nucleatlon  function  Is 
specified  for  the  second  layer.  The  total  number  of  pits  nucle¬ 
ated  at  time  x  is 


n(T) 


T 

6n(T .♦)d^ 

0 


(7) 
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Then  the  volume  removed  from  the  second  layer  at  time  t  due  to 
all  the  pits  present  at  time  t  will  be 


fiVaCtjX)  ■  DO|(t-T)*n(T)  (8) 

where  the  linear  growth  rate  on  the  second  layer  Is  taken  to  be 
different  from  that  on  the  first.  Combining  Eqs.  (6)  to  (8)  and 
Integrating  with  respect  to  x,  we  find  that  the  total  volume 
eroded  from  the  second  layer  at  time  t  Is  simply 


AVa(t) 


0| 


X 

Ia(x-<»)* 


0  0 


(9) 


(dAV,(^|») 

7^ — 


Or  finally,  Eq.  (9)  can  be  put  Into  a  more  convenient  form  by 
differentiation  with  respect  to  t. 


dAVi(t) 

“"IE - 


ft 


2  a| 


(t-x) 


X 

Ij(x-l|>)' 


0  0 


(10) 


^dAV|(<») 

— 


dAV2(<»)^ 

— 3i; — J 


-13- 


Eq.  (10)  Is  slightly  more  complicated  than  Eq.  (4);  but  once 
AViCip)  Is  found  from  Eq.  (4),  It  can  be  inserted  in  Eq.  (10)  for 
evaluating  • 


The  derivation  of  the  equations  governing  material  removal 
from  subsequent  layers  proceeds  in  a  manner  analogous  to  that  for 
the  second  layer.  In  summary,  the  calculation  of  the  weight  loss 
(or  rate  of  erosion)  experienced  by  a  material  body  subjected  to 
an  erosive  environment  requires  solving  the  following  system  of 
equations . 


AVi(t)  -  G?  (t-T)*Ii(T)[V-AVi(T)]dT 


dAV^(t) 


2  G*  (t-T)  Ijjj(t-iI^)* 


\  3^  dhji  y 


dij^dT 


(11) 


where  o  ■  2,3f..«.»N,  and 


AW(t)  -  pAV(t)  -  p  J  AV  (t) 


(12) 


is  the  total  weight  loss  at  any  time  t. 
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The  system  of  equations  has  been  programmed  for  digital 
computer  computations.  These  computations  are  extremely  effi¬ 
cient,  since  a  considerable  simplification  in  their  numerical 
evaluation  was  obtained  through  an  observation  pertaining  to  the 
general  form  of  the  equivalent  system  of  finite-difference 
equations.  There  is  a  considerable  degree  of  flexibility  in  the 
functional  form  of  the  pit  nucleation  probabilities,  the  growth 
rates,  and  mean  pit  depths,  since  the  computer  code  is  not 
restricted  in  this  regard.  The  explicit  details  of  the  devia¬ 
tion  of  Eq.  (11)  can  be  readily  modified  when  experimental 
observations  indicate  that  such  modifications  are  in  order.  The 
computer  program  is  provided  in  the  Appendix. 

n  p\ 

Heymann''  '  developed  a  conceptually- similar  model  which 
incorporated  the  nucleation  and  growth  of  eroding  areas  in  terms 
of  a  log-Normal  frequency  distribution  function.  Heymann  demon¬ 
strated  that  qualitative  agreement  can  be  obtained  with  various 
forms  of  experimental  erosion  rate- time  curves  through  proper 
selection  of  the  parameters  in  the  distribution  function.  His 
model  can  only  describe  the  erosion  data  once  it  is  obtained. 

The  model  developed  here  could  be  used  for  a  similar  purpose. 

The  problem  of  analytically  determining  the  erosion  of  a 
material  which  behaves  according  to  our  model  is  reduced  to  the 
problem  of  specifying  nucleation  and  growth  kinetics  for  erosion 
pits.  The  value  of  this  model's  calculations  as  a  predictive 


tool  In  the  development  of  erosion  resistant  materials  depends 
on  how  much  Information  pertaining  to  the  kinematics  of  the 
erosive  environment  and  the  material  characteristics  of  the  eroding 
solid  can  be  Incorporated  Into  the  explicit  analytic  representa¬ 
tion  of  the  pit  nucleatlon  and  growth  rate  functions.  It  would 
be  satisfying  to  prescribe  these  functions  on  a  purely  analytical 
basis;  however,  this  Is  not  possible  at  the  present  time.  This 
research  considers  the  relevant  analytical  approaches  in  conjunc¬ 
tion  with  a  series  of  controlled  erosion  tests. 
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SOLID  PARTICLE  IMPACTS  ON  BRITTLE  MATERIALS 

i'  ‘  I 

I 

4  t 

Itie  statistical  parameters  required  In  ^he  erosion  theory 

I  . 

developed  In  the  previous  section  can  ‘be  obtained  by  the  direct 

I  1  ■  • 

observation  of  pit  nucleatlon  and  growth  In  materials  which 
erode  In  this  manner.  However,  the  predictive 'character  of  the 
model  Is  greatly  enhanced  If  the  kinematic  parameters  associated 

>  I  !  '  ' 

with  specific  erosion  conditions  and  material  properties  for  general 
material  classes  can  be  Introduced  Into  the  general  formulation. 

In  order  to  accomplish  this  objective,  the  general , erosion  model  ' 

» 

Is  applied  to  the  simplest  conditions  of  solid  particle  Impacts 
on  brittle  materials.  The  Initial  evaluation  of -the  erosion  model. 

I  • 

will  be  balsed  on  the  following  assumptions: 

(  •  I 

It.  ,  I  ■  , 

(1)  The  eroded  solid  can  be  represented  by  a  semi- Infinite 

half' space  which  Initially  has  a  perfectly  plane  surface^  ' 

(2)  The  eroded  solid  Is  a  linearly  elastic  material' 
which  Is  homogeneous  and  Isotropic. 

*  I  1  I  •  ■  I  I 

(3)  The  eroding^  particles  are  spherical  In  form. 

(4,)  The  eroding  parvlcles  are  homogeneous,  and  Isotropic 

*  I 

linearly  elastic  solids.  , 

I 

(5)  The  Impact  velocities  are  low  enough  that  the  Hertzian 

I  .  ■ 

theory  of  Impact  is  applicable.  '  , 

I 

(6)  The  eroding  particles  only  collide,  wll^h  the  eroded 
solid  In  a  direction  which  Is  perpendicular  ,to  Its  Initially 

■  I  • 

plane  surface. 
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In  order  to  make  the  analysis  of  the  particle  Impacts  as 
manageable  as  possible.  It  Is  desirable  to  work  within  the  domain 
of  the  Hertzian  theory  of  Impact  which  Is  one  of  the  few  analytical 
approaches  to  Impact  phenomena  that  can  be  handled  without  the  aid 
of  elaborate  computer  codes.  This  choice  allows  us  to  take 
advantage  of  the  extensive  prior  theoretical  and  experimental 
work  on  Hertzian  Impacts  of  purely  elastic  bodies. 

Inorganic  glasses  were  considered  to  be  most  representative 
of  the  Idealized  brittle,  elastic  materials  required  In  the 
Hertzian  analysis.  For  the  experimental  studies  of  the  erosion 
mechanisms  and  erosion  characteristics  under  the  prescribed 
conditions,  fused  silica  (Corning  79^0)  was  selected  for  the 
glass  specimens  and  soda  lime  glass  beads  are  being  used  for  the 
Impacting  particles.  The  criteria  for  fracture  Initiation  of 
glass  surfaces  loaded  by  spherical  Indenters  Is  also  considered 
In  this  section. 

A.  HERTZIAN  THEORY  OF  IMPACT 

Since  much  of  the  discussion  which  follows  is  based  on 
the  results  of  the  Hertizan  theory  of  Impact,  the  pertinent 
aspects  of  this  theory  will  be  presented  for  future  reference. 

Hertz assumed  that  the  collision  of  two  bodies  could  be 
regarded  as  a  statical  problem;  the  compression  at  the  place  of 
contact  Is  regarded  at  gradually  produced  and  as  subsiding  completely 
by  reversal  of  the  process  by  which  It  was  produced.  Hertz  theory 


-18- 


of  Impact  for  perfectly  elastic  bodies  does  not  account  for 
dissipation  of  energy;  therefore,  no  loss  In  the  total  kinetic 
energy  of  the  two  bodies  will  occur.  This  restriction  has  been 
relaxed  In  more  recent  analyses  of  the  Impact  of  a  spherical 
Indenter  on  a  viscoelastic  half-space. 

The  problem  originally  considered  by  Hertz  was  to 
determine  for  the  case  of  two  elastic  Isotropic  bodies  which 
touch  each  other  over  a  very  small  part  of  their  surface  and 
exert  upon  each  other  a  finite  pressure,  distributed  over  the 
common  area  of  contact,  the  form  of  the  common  surface,  the  dis¬ 
tribution  of  pressure  over  It,  and  a  relation  between  the  distance 
the  bodies  approach  each  other  and  the  pressure  existing  on  the 
contacting  surface  of  each  body.  The  results  given  here  will  be 
restricted  to  the  case  of  a  spherical  elastic  body  Indenting  an 
elastic  half-space.  A  detailed  account  of  the  stress  distribution 
for  this  case  was  first  given  by  Huber the  specialized  theory 
Is  often  referred  to  as  the  Hertz-Huber  theory. 

Hertz  makes  the  following  simplifying  assumptions: 

(a)  The  surface  of  contact  Is  always  small  In 
comparison  with  the  dimensions  of  the  bodies  Involved. 


(b)  The  surfaces  in  contact  are  imagined  as  perfectly 
smooth,  i.e.  we  assume  that  only  a  normal  pressure  acts  between 
the  parts  in  contact. 

To  apply  the  static  theory  to  the  impact  of  perfectly  elastic 
bodies  the  following  additional  assumptions  were  made  by  Hertz: 

(c)  ...the  elastic  state  of  the  two  bodies  near  the 
point  of  Impact  during  the  whole  duration  of  impact  is  very  nearly 
the  same  as  the  state  of  equilibriiun  which  would  be  produced  by 
the  total  pressure  subsisting  at  any  Instant  between  the  two 
bodies,  supposing  it  to  act  for  a  long  time. 

Thus,  no  account  of  vibrations  present  in  the  bodies  due  to  impact 
was  undertaken,  and  the  pressure  between  the  two  bodies  is  assumed 
to  be  gradually  applied  even  though  the  entire  time  the  bodies  are 
in  contact  is  relatively  short. 

(d)  ...the  time  of  impact  is  also  large  compared 
with  the  time  taken  by  elastic  waves  to  traverse  the  impinging 
bodies  from  end  to  end. 

The  problem  under  consideration  becomes  that  of  an 
clastic  body  subjected  to  a  pressure  distributed  over  a  very  small 
portion  of  its  surface.  Hertz  thus  reduced  the  problem  to  the 
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the  case  of  a  body  bounded  by  an  infinite  plane  which  is  appropriate 
if  the  dimensions  of  the  body  are  lar<;e  compared  with  those  of  the 
area  of  the  load. 


The  radius  a  of  the  circle  of  contact  and  the  distance 


of  approach  o 

of  the 

two  bodies  are  then  given  by 

1/3 

a 

- 

PR  (k,^  +  k^)) 

(13) 

a 

(14) 

where  P  is  the  normal  load  on  the  indenter,  R  is  the  Indenter 
radius,  and 


Kl  .  Iz^ll  .  Kj  .  (15) 

are  the  elastic  moduli  for  the  sphere  and  half-spwe,  respectively, 
(  V  is  Poisson's  ratio  and  E  is  Young's  modulus).  If  the  magnitude 
of  the  maximum  pressure  is  scaled  to  the  radius  of  the  contract 
area,  the  pressure  distribution  over  the  contact  area  is  shown  to 
be  hemispherical  for  the  particular  case  under  consideration. 

Having  the  magnitude  of  the  contact  surface  and  the 
pressure  acting  on  it,  the  elastic  stress  distribution  in  the 
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half- space  can  be  evaluated  Simple  formulas  for  the  stress 
components  can  only  be  obtained  along  the  axis  of  symmetry  (the 
positive  s-axls  directed  Into  the  Interior  of  the  half-space)  and 
over  the  plane  surface  of  the  half-space.  The  determination  of 
the  stress  components  at  other  points  In  the  half-space  requires 
numerical  analysis.  The  radial  tensile  stress  on  the  surface  of 
the  half-space  outside  the  contact  area  Is 


Or 


P 


(r-  a) 


(16) 


The  radial  stress  Is  compressive  for  r  -  a.  TTie  circumferential 
component  of  the  stress  os  Is  compressive  over  the  entire  surface 
and  of  the  same  magnitude  as  for  r  -  a.  Along  the  z-axls 
(axis  of  symmetry)  °  x  compressive  between  the 

surface  and  a  critical  depth  which  depends  on  v  .  At  this  critical 
depth  both  o^and  vanish  and  are  tensile  for  all  greater  depvns. 
The  0^  stress  component  Is  compressive  at  all  depths  along  the 
axis. 


The  maximum  radial  tensile  stress  Is  located  at  the 
circular  boundary  of  the  contact  surface,  Bq.  (16),  and  has  a 
magnitude  of 


'max 


1-2  Vg 


fl/3  r-2/3 


(r  -  »)(17) 
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To  study  the  Impact  of  a  spherical  body  on  an  elastic 


surface,  the  previous  relations  from  the  Hertzian  theory  are 
restated  In  terms  of  the  velocity  of  approach  v  Instead  of  the 
load  P  which  Is  now  an  unknown.  In  order  to  establish  a  relation 
between  P  and  v,  the  kinetic  energy  of  the  Impacting  sphere  Is 
equated  to  the  strain  energy  at  the  time  of  maximum  Impression. 


wnere 
Eq.  (U), 


1/^  3x  2  . 

p,it  R  '  V 


max 


P(a)aw 


is  the  density  of  the  Impacting  sphere. 


Using 


The  duration  of  contact  T  between  the  two  Impacting  bodies  Is 
given  by 


T 


2.9^3  C  I 


2/5 

TiPj  (kj  +  kg)] 


(19) 


The  results  of  the  Hertzian  theory  of  Impact  have 
been  Investigated  experimentally  since  the  time  of  Hertz  and  have 
been  confirmed  for  a  broad  range  of  materials  and  Impact  conditions 
The  greatest  deviations  from  the  Hertlzan  predictions  occur  In  the 
case  of  very  small  diameter  spheres  Impacting  on  a  plane  surface 
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which  l8  p*rticularly  relevant  to  solid  particle  erosion.  As  the 
duration  of  contact  decreases*  either  due  to  Impact  by  small 
particles  or  high  Impact  velocities,  the  dynamic  effects  of  the 
lmp<'\ct  must  be  taken  Into  account. 


B.  FRACTURE  OF  GLASS  SURFACES  BY  SPHERICAL  INDENTERS 

In  a  theory  of  erosion  for  brittle  materials  through 
repeated  Hertzian  Impacts,  the  dominant  consideration  is  how  the 
Hertzian  stress  distribution  influences  the  nucleation  and  pro¬ 
pagation  of  brittle  fractures.  If  an  elastic  sphere  is  pressed 
onto  en  elastic  surface  containing  a  dense  distribution  of  fine 
flaws  under  steadily  increasing  load,  no  appreciable  crack  propa¬ 
gation  is  observed  until  the  load  exceeds  a  critical  magnitude  P^. 
The  evaluation  of  the  critical  fracture  load  for  a  spherical 
indenter  pressed  against  a  glass  plate  has  been  an  active  area  of 
research  since  the  time  of  Hertz. 

The  crack  produced  quasistatically  in  a  large  glass 
plate  by  a  spherical  indenter  appears  as  a  ring  on  the  surface 
and  spread  out  below  the  eurface  to  form  the  fmstrum  of  a  cone. 
Once  a  fracture  is  formed,  however,  the  stress  distribution  is 
modified,  since  the  fracture  surfaces  act  as  free  boundaries  for 
which  the  normal  stress  is  zero:  the  Hertz-Huber  analysis  will  no 
longer  be  valid.  In  this  regard,  it  is  noted  however  that  Sheldon 
and  Finnie^^^  have  provided  experimental  evidence  that  the  Hertzian 
expression  for  the  distance  of  approach  betw^^en  a  1/6  in.  diameter 
steel  sphere  impressing  a  pyrex  glass  plate  is  not  appreciably  in 
error  even  after  the  application  of  loads  in  excess  of  120  times 
the  initial  fracture  load  and  multiple  ring  fractures  were  formed. 
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The  Hertz-Huber  equations  give  an  accurate  representation 

of  the  stress  distribution  In  glass  plates  before  fracture 

Initiates.  Prom  Eq.  (17)  the  maximum  radial  tensile  stress  Is 

located  at  the  circular  boundary  of  the  contact  surface  for  the 

Inhomogeneous  stress  distribution.  In  the  case  of  brittle 

materials,  such  as  glass,  a  simple  fracture  criterion  Is  that 

failure  will  occur  when  a  critical  value  of  the  maximum  tensile 

stress  Is  reached.  If  the  tensile  strength  of  glass  Is  assumed 

to  be  constant,  then,  according  tc  Eq.  (17),  the  load  P.  required 

to  produce  cracking  should  vary  with  the  square  of  the  Indenter 

2 

radius.  This  relation,  P.  proportional  to  R  ,  has  been  confirmed 
experimentally  for  large  sphere  radii  (minimum  radii  on  the  order 
of  one  Inch).  However,  Auerbach^^^  found  that  as  the  radius  of 
the  Indenter  decreased,  the  magnitude  of  the  critical  load  for 
fracture  Increased  In  direct  proportion  to  R.  The  relation 

P-^  -  H  (a  constant)  (20) 

has  been  designated  Auerbach's  law  by  Roesler^^^^  Substitution 
of  this  relation  Into  Eq.  (17)  means  that  the  fracture  stress 
varies  as  the  one-third  power  of  the  Indenter  radius.  The 
observed  size  effect  In  the  fracture  of  glass  has  led  to  consider¬ 
able  controversy  as  to  the  physical  explanation  of  Auerbach's  luw. 


-26- 


One  explanation  of  thla  size  effect  In  glass  Is 
based  on  modification  of  Griffith's  criterion  for  brittle  fracture 
applied  to  flaws  In  glass  surfaces.  The  modification  takes  Into 
account  the  fact  that  the  stress  required  to  Initiate  a  fracture 
from  a  particular  flaw  In  brittle  materials  depends  on  the  size, 
location,  and  orientation  of  the  flaw  In  the  Inhomogeneous  stress 
field  produced  by  a  spherical  Indenter  on  a  plate.  It  Is  then 
conjectured  that  a  larger  stressed  area  (larger  Indenter  radius) 
has  a  greater  probability  of  Including  a  particularly  weak  flaw. 
However,  as  Roesler^^^^  points  out,  according  to  this  reasoning 
the  scatter  In  the  evaluation  of  the  critical  load  should  Increase 
as  the  size  of  the  radius  of  the  Indenter  decreases.  This  scatter 
Is  not  observed  experimentally;  the  scatter  In  the  evaluation  of 
the  critical  load  Is  the  same  for  small  as  well  as  large  Indenters. 
Statistical  theories  of  fracture  from  assumed  surface  flaw  distri¬ 
butions  In  a  Hertzian  stress  field  have  been  advanced  recently  by 
Tsai  and  Kolaky^^^^  Oh  and  Flnnle^^^\  and  Fisher^ Fisher 
demonstrates  that  a  particular  form  of  flaw  distribution  Is  required 
to  obtain  the  Auerbach  relation,  Eq.  (20).  It  Is  highly  Improb¬ 
able  that  special  forms  of  the  flaw  distribution  are  present  In 
all  of  the  glass  samples  tested. 

An  alternative  theoretical  treatment  of  the  problem 

( li|^ 

has  been  developed  by  Frank  and  Lawn^  ,  who  employ  Roesler's 

(10) 

energy  scaling  criterion  .  Roesler  showed  that  the  various 


-27- 


empirical  relations  for  the  production  of  cone  cracks  can  be  put 
into  the  form  that  Just  before  fracture  the  elastically  stored 
energy  is  proportional  to  the  area  of  the  circle  of  contact. 

Since  cone  cracks  are  geometrically  similar,  scaling  to  the  circle 
of  contact,  Roesler  suggested  that  .he  fundamental  form  of 
Auerbach’s  law  is  the  proportionality  between  the  elastically 
stored  energy  and  the  crack  surface  area.  Frank  and  Lawn^^^^ 
studied  conditions  for  the  Instability  of  embedded  cracks  of  the 
Griffith  type  in  the  inhomogeneous  Hertzian  stress  field.  They 
developed  a  procedure  for  determining  stress  intensity  factors 
at  the  tips  of  these  flaws  and  the  energy  release  on  incremental 
crack  propagation.  These  relations  were  used  to  obtain  criteria 
for  Hertzian  fracture  which  predict  behavior  substantially  in 
accord  with  experiment  and  lead  to  estimates  of  the  surface  tension 
of  fracture  surfaces  in  glass  which  are  in  reasonable  agreement 
with  those  obtained  through  other  techniques. 

Lemgltan  and  Lawn^^^^  cranpared  the  predictions  of  the 
flaw  statistical  theories  and  the  energy  balance  theory  of  Frank 
and  Lawn  for  abraded  glass  surfaces  which  can  be  assumed  to  contain 
a  uniform  distribution  of  flaws.  They  demonstrate  that  the 
experimental  evidence  does  not  conform  to  the  results  of  the  flaw 
statistical  theory  for  this  particular  case  which  provides  an 
explicit  representation  for  the  flaw  distribution.  According  to 
the  flaw  statistical  theories  there  should  be  no  difficulty  in 
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encountering  a  critical  flaw;  however,  it  is  found  experimentally 
that  the  fracture  load  is  greater  than  that  evaluated  for  as- 
received  glass  surfaces. 

Results  from  the  energy  balance  approach  are  being 
considered  in  conjunction  with  establishing  fracture  criteria  for 
the  Impact  conditions  which  are  of  Interest  in  erosion  processes. 
The  production  of  ring  fractures  on  the  glass  surfaces  by  multiple 
particle  impacts  is  the  elemental  fracture  mechanism  in  this 
study,  however  a  single  ring  fracture  may  or  may  not  be  the 
nucleus  for  a  growing  erosion  pit.  The  nucleatlon  of  an  erosion 
pit  may  occur  through  one  of  several  processes  whose  relative 
probabilities  are  difficult  to  assess  without  evaluating  a  detailed 
model.  Thus,  the  critical  step  in  the  nucleatlon  of  an  erosion 
pit  may  be  the  impact  of  a  sphere  causing  a  ring  fracture  in  the 
immediate  neighborhood  of  a  previously  formed  ring  fracture,  it 
may  require  three  adjacent  ring  fractures,  and  more  complex 
nucleating  steps  may  easily  be  Imagined.  Thus  the  probability  of 
pit  nucleatlon  as  a  function  of  time  of  exposure  to  erosion  must 
be  determined  through  an  experimental  study  of  the  dominant  pit 
nucleatlon  mechanism  in  conjunction  with  the  analyses  of  Hertzian 
Impact.  The  experimental  study  is  described  in  the  next  section. 
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V.  DESCRIPTION  OF  EROSION  EXPERIMENTS 


The  experimental  program  Is  designed  to  evaluate  the 
damage  to  the  eroded  surface  of  fused  silica  specimens  after 
Incremental  exposures  to  the  direct  Impact  of  spherical  glass 
beads  In  the  AFML-Bell  erosion  apparatus.  A  procedure  Is 
established  for  evaluating  the  statistical  parameters  required 
In  the  theory,  and  the  fundamental  mechanisms  for  the  nucleatlon 
of  pits  In  fused  silica  are  determined  experimentally. 

A.  EXPERIMENTAL  PROCEDURES 

Due  to  the  low  Impact  velocities  used  In  this  study, 

It  was  possible  to  design  specimens  with  a  sizable  flat  surface 
exposed  to  the  eroding  particles.  This  geometry  eliminated  the 
effect  of  the  variable  angles  of  attack  associated  with  the  aero¬ 
dynamic  shapes  usually  employed  In  the  AFML-Bell  erosion  apparatus. 

A  special  sample  holder  was  fabricated  for  the  rotating  arm  so 
that  these  flat-faced  specimens  could  be  used.  Fused  silica 
(Corning  79^0)  was  selected  for  the  erosion  targets.  The  general 
specimen  dimensions  are  Indicated  In  Figure  2.  The  diameters  of 
the  glass  beads  (Impacting  particles)  were  selected  so  that  they 
would  be  small  In  comparison  with  lateral  dimensions  of  the  specimen. 
Under  these  conditions  the  effects  of  the  Impacts  will  be  quite 
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localized;  the  Influence  of  the  finite  boundaries  of  the  specimen 
will  be  minimized,  and  the  central  region  of  the  specimen  would 
approximate  a  mathematical  half-space  as  required  In  the  Hertzian 
analysis.  In  order  to  provide  an  accurate  appraisal  of  the  model, 
spherical  glass  beads  (3M  Superbrlte-Knoop  Hardness  No.  525)  with 
diameters  of  2.8  mils  (70^  )  and  11.4  mils  (290^)  were  used.  The 
beads  showed  some  variation  In  size  from  these  dimensions,  but  the 
size  distribution  was  quite  narrow  around  the  diameter  Indicated. 

A  comparison  of  the  two  bead  sizes  and  the  extent  of  their  spheri¬ 
city  can  be  seen  In  Figure  3*  The  mechanical  properties  of  the 
fused  silica  (Coming  7940)  and  the  glass  beads  are  provided  In 
Table  I. 

The  test  conditions  for  the  fused  silica  samples 
are  given  In  Table  II.  Particl''  Impact  velocities  of  200  and 
300  fps  were  selected  on  the  basis  of  adjunct  experimentation  using 
pyrex  samples.  Impact  velocities  were  required  which  would  be 
capable  of  producing  ring  fractures  on  the  surface  of  the  fused 
silica  but  were  not  of  a  magnitude  that  the  Initial  and  Inter¬ 
mediate  stages  of  the  erosion  process  would  be  completed  before 
the  required  statistical  data  were  obtained. 

In  the  erosion  experiments  conducted  In  the  AFML-Bell 
erosion  apparatus  shown  In  Figure  4,  the  glass  beads  are  taken 
from  the  sand  hopper  by  a  screw-type  feeder  and  fall  down  the 
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TABLE  I 

MECHANICAL  PROPERTIES  OF 
FUSED  SILICA  AND  SODA  LIME  GLASS 


Soda  Lime  Glass 


Density  (Ib./ln.^) 
Poisson's  Ratio 
Young's  Modulus  (Ib/ln^) 


0.0866 

0.22 

10.0  X  10^ 


Fused  Silica 
0.0795 

0.17 

10.68  X  10^ 
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TABLE  II 


EROSION  TEST  SCHEDULE 


SPECIMEN 

BEAD  SIZE 
(MILS) 

BEAD  FLOW 
RATE 

^LB/MIN. ) 

IMPACT 
VELOCITY 
(FT. /SEC) 

G1 

2.8 

2 

200 

02 

11.4 

1 

200 

03 

11.4 

1 

300 

-35- 


-36- 


1  3/8  Irch  diameter  blast  tube  which  Is  approximately  7  feet  long. 
They  are  dropped  through  the  plane  of  the  rotating  arm  from  an 
elliptical  nozzle  which  has  a  major  axis  of  1  7/8  Inch  and  a 
minor  axis  of  7/8  Inch.  Due  to  the  elliptical  opening  at  the  end 
of  the  drop  tube  the  distribution  of  beads  Is  nonuniform  along 
the  length  of  the  specimen.  Modification  of  the  nozzle  Is  under 
consideration  to  provide  a  uniform  distribution  of  particles  along 
the  length  of  the  specimen. 

The  rate  of  flow  of  the  beads  was  calibrated  by  weigh¬ 
ing  the  beads  which  fall  from  the  tube  for  a  number  of  settings  on 
the  feed  mechanism.  Although  It  was  not  possible  to  ascertain 
directly  the  actual  number  of  beads  striking  the  surface  of  the 
specimen,  a  theoretical  calculation  of  the  maximum  number  of 
Impacts  Is  possible. 

Let  N  be  the  number  of  beads  released  from  the  hopper 
per  second  and  x  be  the  distance  a  bead  falls  from  the  feeder 
due  to  gravity  at  time  t, 

X  .  (21) 

The  density  of  the  stream  of  beads  will  vary  along 
the  length  of  the  tube.  At  a  distance  x  from  the  release  point, 
the  spatially- dependent  denslty^(x)  Is  given  by 

./..X  llm  N(t4At)  -  Nt  11m  N/Axx-1 

'l<  a - 

or 
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(x) 


N 


(22) 


using  Eq.  (21).  A  Is  the  cross-sectional  area  of  the  stream  of 
beads.  The  maximum  number  of  impacts  on  the  specimen  per 
revolution  of  the  rotating  arm  is 

beads  \  /volume  intercepted^ 


N, 


/  pecQB  \  /voiuine  inuercepteax 

'unit  volume'  '  revolution  ' 


"  h 


JL_)  (Ad)  .  ^ 


(23) 


where  x  ■  h  is  the  distance  frora  the  feeder  to  the  plane  of  the 
rotating  arm  and  d  is  the  height  of  the  specimens. 


The  number  of  11  mil  beads  In  1  lb.  was  determined  to 
7  Q  ^ 

be  1.5  X  10  beads.  There  are  10^  2.8  mil  beads  in  1  lb.  Thus 
for  typical  dimensions  of  h  «  7  ft.,  d  ■  0.0625  ft.,  and  a  flow 
rate  of  1  Ib/mln.,  Nj  ■  ^9»000  impacts/rev.  (2.8  mil  beads)  and 
Nj  “7^0  Impacts/rev.  (11  mil  beads). 

The  velocity,  v,  of  the  specimen  near  its  midpoint 
on  the  rotating  arm  is  related  to  the  rotational  speed  of  the 
arm,i7,  measured  in  revolutions  per  minute  (RPM)  by 

n  -  i.o€v 


when  V  is  measured  in  ft/sec. 
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The  basic  premlae  of  the  theory  developed  In  Section 
III  la  that  erosion  occurs  through  preferential  growth  of  eroding 
regions  on  the  surface  of  the  specimen.  The  pit  nucleatlon 
function  and  mean  growth  rate  are  the  statistical  parameters 
which  represent  this  process.  A  straight-forward  procedure  for 
obtaining  explicit  representations  of  these  functions  is  the 
direct  observation  of  nucleated  erosion  sites  on  the  sample 
surface  as  a  function  of  time.  The  greater  the  number  of  sites 
tabulated,  the  more  accurate  the  statistical  data  will  be.  The 
experimental  program  Is  designed  to  evaluate  the  damage  to  the 
eroded  surface  of  fused  silica  specimens  after  Incremental 
exposures  to  the  direct  impact  of  spherical  glass  beads  in  the 
AFML-Bell  erosion  apparatus. 

The  specimens  were  handled  with  cotton  gloves,  weighed 
and  carefully  mounted  in  the  specimen  holder  on  the  blade.  They 
were  exposed  to  the  direct  Impacts  of  the  glass  beads  for  time 
periods  which  increased  with  the  increase  in  the  cummulatlve 
time  of  exposure.  The  damage  process  in  each  specimen  was  monitored 
microscopically  after  each  time  increment  in  order  to  evaluate 
the  statistical  parameters  for  the  proposed  theory  and  determine 
the  erosion  mechanisms  responsible  for  the  nucleation  of  pits 
In  the  eroding  surface.  The  scale  of  the  damage  was  such  that 
optical  microscopy  could  be  used  to  make  the  necessary  observa¬ 
tions.  A  special  staging  plate  was  designed  and  fabricated  for 
the  Bausch  and  Lomb  Metallograi^  which  made  it  possible  to  return 
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conslatently  to  within  a  few  mils  of  a  particular  location  on  the 
surface  of  a  specimen.  A  magnification  of  IlOX  was  found  to  be 
adequate  for  observing  the  details  of  the  erosion  process  which 
were  seen  best  using  dark  field. 

The  experimental  procedure  for  recording  the  extent 
of  the  erosion  damage  after  each  Incremental  exposure  on  the 
rotating  arm  was  first  to  survey  the  entire  surface  for  any 
unusual  feature  produced  by  the  particle  Impacts.  Then  a  total 
of  seventy-five  sites,  approximately  sixty  nils  In  diameter 
located  In  the  central  region  of  the  specimen  are  Investigated 
In  detail.  In  the  Initial  stages  of  the  erosion  process  the 
number  of  ring  fractures  Is  counted  for  each  site  along  with  the 
number  of  double,  triple,  and  more  complex  ring  fracture  Inter¬ 
actions.  When  possible,  the  observed  damage  was  recorded  photo¬ 
graphically.  The  depth  of  a  particular  feature  on  the  surface  can 
be  measured  to  within  a  few  microns  using  the  Metallograph.  By 
setting  the  focus  at  a  level  below  the  surface.  It  Is  possible 
to  note  the  extent  of  the  damage  at  this  level  based  on  areas 
which  are  In  focus  as  the  surface  of  the  specimen  Is  scanned. 

It  Is,  therefore,  possible  to  divide  the  specimen  Into  a  number 
of  layers  using  the  optical  microscope  In  conformity  with  the 
proposed  model. 
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In  order  to  minimise  the  effect  of  surface  imperfections, 
three  specimens  of  fused  silica  (01  to  03)  were  prepared  with 
highly-polished  surfaces.  Replicas  were  made  of  the  surfaces  and 
viewed  in  the  transmission  electron  microscope  at  magnifications 
up  to  38,000X.  The  highly- polished  surfaces  were  featureless: 
no  surface  flaws  could  be  detected  and  only  rarely  was  an  inclusion 
observed.  As  indicated  in  the  literature  surface  flaws 

cannot  be  observed  in  this  way.  Other  means  devised  for  the 
direct  observations  of  flaw  distributions  in  glass  surfaces  are 
still  viewed  with  uncertainty  as  to  the  evaluation  ‘^of  the  actual 
distribution  and  that  which  was  induced  by  the  experimental  pro¬ 
cedure  used.  In  this  initial  study  no  detailed  investigations  were 
made  of  the  flaw  distributions  on  the  original  surfaces  of  the 
specimens,  since  the  perfection  of  the  highlj'^-polished  surfeu^es 
would  be  disrupted  by  the  experimental  procedures.  When  flaw 
distributions  are  required,  test  bars  of  fused  silica  which 
received  the  same  treatment  as  the  specimens  are  available  for 
this  evaluatim  either  through  direct  or  Indirect  means. 

The  equivalent  static  load,  radius  of  contact  area, 
and  duration  of  iiapact  can  be  calculated  from  the  Herts  theory, 

Eq.  (l8),  (13)  and  (19 )»  using  the  velocities  for  the  erosion  test 
conditions  specified  in  Table  II  and  the  mechanical  properties 
in  Table  I.  The  results  of  these  computations  are  listed  in 
Table  III. 
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TABLE  III 


HERTIZAN  IMPACT  PARAMETERS  FOR 
FUSED  SILICA 


2.8  Mil 
Glass 
Beads 

Equivalent  Static 
Load,  P(lb.) 

Radius  of  Contact 
Circle,  a(ln.) 

Duration  of 
Contact,  T(sec) 

200  fps 

0.215 

0.350  X  10"3 

1.19  X  10" 

300  fps 

0.355 

0.415  X  10"3 

1.10  X  10" 

11  Mil 
Glass 
Beads 


200  fps 

3.56 

1.430  X  10"3 

4.76  X  10' 

300  fps 

5.88 

1.692  X  10"3 

4.40  X  10' 
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The  results  In  Table  III  show  that  the  maximum  contact 


area  Is  large  In  comparison  to  the  diameter  of  the  Impacting  beads 
which  cast  some  doubt  on  the  satisfaction  of  the  basic  assumption 
of  the  Hertzian  analysis.  Assuming  the  Hertzian  approach  is 
valid,  the  calculated  duration  of  contact  (on  the  order  of  10  ' 
seconds)  indicates  that  a  dllatational  wave*  will  travel  a  maxi¬ 
mum  on  0.1  In.  from  the  initial  point  of  contact  during  this 
time  Interval.  The  shear  and  Rayleigh  waves  will  have  traveled 
only  0.06  Inch.  On  the  basis  of  these  distances  In  relation  to 
the  dimensions  of  the  maximum  contact  area  and  bead*"  size,  the 
quaslstatlc  loading  requirements  for  the  treatment  of  Impact 
phenomena  can  also  be  questioned.  The  Impact  conditions  in  the 
experimental  program  were  dictated  by  the  size  of  the  sample 
which  could  be  accommodated  on  the  rotating  arm.  This  placed  a 
restriction  on  the  size  of  the  impacting  particles. 


*The  speed  of  a  dllatational  wave  In  fused  silica  is  19,600  ft/sec 
at  70®P“  and  a  shear  wave  travels  at  12,350  ft /sec. 
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B.  EROSION  MECHANISMS  IN  FUSED  SILICA 


Although  a  variety  of  inorganic  glasses  have  been 
used  in  erosion  studies  reported  in  the  literature,  no  information 
could  be  found  on  the  operative  erosion  mechanisms  associated  with 
the  various  stages  of  the  weight-loss  data  for  glasses  subjected 
to  multiple  solid  particle  impacts.  The  erosion  data  for  solid 
particles  impacting  brittle  materials  are  generally  reported  in 
terms  of  plots  of  the  mass  of  specimen  loss  per  unit  mass  of 
eroding  particles  (or  an  equivalent  measure  of  erosion)  as  a 
function  of  the  particle  impingement  angles.  Such  an  approach  only 
provides  the  steady- state  erosion  rate  for  a  prescribed  particle 
size  (or  distribution  of  particle  sizes),  particle  shape,  and 
impact  velocity.  The  work  reported  here  investigates  the  erosion 
mechanisms  in  an  inorganic  glass  when  eroded  by  multiple  collisions 
with  relatively  small  particles. 

Specimen  G1  was  exposed  to  impacts  by  2.8  mil  beads  at 
200  fps  (or  the  equivalent  of  212  RPM  on  the  rotating  arm).  The 
bead  flow  rate  was  2  Ib/min.,  so  according  to  Eq.  (23)  the  maxi¬ 
mum  number  of  bead  impacts  should  be  on  the  order  of  98,000  impabts/ 
rev.  The  specimen  had  a  cummulative  exposure  time  of  72  minutes. 

The  erosion  test  for  specimen  01  was  exploratory. 

During  the  early  stages  of  erosion,  under  8  min.,  the  eroded 
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surface  was  composed  of  a  distribution  of  ring  fractures,  small 

»  '  '  I  ' 

pitted  areas  less  than  1  mil  In  diameter,  and  a  distribution  of 
\  \  •  * 

very  small  Imprints,  less  than  0.1  mil,  due  'to  collisions  which 

I 

did  not  produce  measurable  damage.  Tlie  average  diameter  of  the 

I  I  *• 

ring  fractures'  was  3  mils  with  very  little  deviation  from  this 
value.  The  formation  of  these  rfnga  on  the  surface  and  their  ■  ' 

Interaction  was  assumed  to  be  the  primary  mechanism  of  material 

removal.  A  portion  ,of  the  ring  fractures  are  accompanied  by 

1 

material  removal  Just  outside  the  ring  perimeter:  the, Interior 
of  the  circular  ring  Is  undisturbed.  Larger-scale  material 

I 

removal  occurs  In  the  Interstitial  areas  of  non- overlapping  d^ouble 
and  triple  ring  fractures.  A  limited  number  of  more  complex  forma 

‘  ^  I 

of  ring  Intersections  are  beginning  to  develop  at  this  time. 

I 

* 

The  volume  removed  from  the  Interstitial  region 
between  ring  fractures  was  considered  to  be  the  nucleatlon  site 

I 

for  a  pit  on  the  surface.  The  location  of  these  sites  In  the 

central  region  of  the  specimen  were  recorded.  The  larger  pits 

.  ' 

nucleated  were  re-examined  after  subsequent  Incremental  exposures 
to  the  eroding  environment,  and  It  was  found  that  none  of  the 

i  * 

previously  established  pits  showed  any  growth.  The  density  of 
3  mil  ring  fractures  on  the  surface  of  the  fused  silica  as  a 
function  of .exposure  time  Is  a  significant  factor  In  the  continual 
growth  of  pitted  regions  as  defined  above.  The  graph  In  Figure  3 

>  I 

shows  that  the  ring  density  (evaluated  for  the  central  region  of 
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the  specimen)  reaches  a  constant  value  after  an  exposure  time  of 
24  minutes.  Tiils  was  a  somewhat  xinexpected  result.  On  the  basis 
of  the  number  of  Impacts  per  revolution*  the  specimen  will  have 

o 

collided  with  1.66  x  10  beads  after  6  min.  and  three  times  this 
number  at  24  minutes.  This  makes  the  probability  for  a  3-nll  ring 
fracture  to  occur  quite  low.  After  an  exposure  time  of  32  min.  the 
microscopic  observations  made  It  quite  clear  that  the  Interaction 
between  the  3-mll  rings  was  not  the  primary  erosion  mechanism* 
but  that  flne> scale  chipping  of  the  surface  was  the  dominant  erosion 
mechanism. 

The  series  of  micrographs  In  Figure  6  Illustrate  the 
time  evolution  of  the  erosion  process  In  the  vicinity  of  a 
grouping  of  3  rail  ring  fractures.  At  8  minutes  a  number  of  ring 
fractures  have  Interacted  to  remove  a  portion  of  the  Interstitial 
material.  The  faint  outline  of  a  number  of  Isolated  ring  fractures 
can  also  be  seen  In  the  general  area  of  this  flaw.  It  was  of 
Interest  to  see  If  subsequent  Impacts  In  the  vicinity  of  this 
nucleated  site  would  extend  the  area  of  the  material  removed  In 
accordance  with  the  proposed  model.  The  micrograph  of  the  same 
region  at  16  minutes  shows  evidence  of  a  large  number  of  particle 
Impacts*  but  essentially  no  new  ring  fractures  result.  Instead, 
the  solid  particle  collisions  are  restrl*:ted  to  fine-scale  chipping 
which  Is  clearly  seen  to  be  preferential  In  the  vicinity  of  the 
existing  ring  fractures.  Some  of  the  ring  fractures  which  are 


(a)  8  MINUTES 


(b )  16  MINUTES 


FIGURE  6  EROSION  IN  THE  VICINITY  OF  3-MIL  RING 
FRACTORES  ON  SPECIMEN  G1  (IlOX). 
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(c)  32  MINUTES 


(d)  56  MINUTES 


FIGURE  6  EROSION  IN  THE  VICINITY  OF  3-MIL  RING 
FRACTURES  ON  SPECIMEN  G1  (110  X). 
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FIGURE  6 


)  72.  MINUTES 


EROSION  IN  THE  VICINITY  OF  3-MIL  RING 
FRACTURES  ON  SPECIMEN  G1  (IlOX). 
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very  faint  In  Figure  6a  are  now  well-decorated.  As  the  exposure 
time  increases  very  small  pieces  are  preferentially  removed  frc»n 
the  regions  of  initial  damage.  The  original  interstitial  region 
between  ring  fractures  exhibits  no  growth  by  subsequent  ring 
fracture  interactions.  Material  is  continually  removed  by  chipping 
of  material  from  the  periphery  and  within  the  eroded  regions, 
rhis  process  continues  until  the  original  surface  of  the  sample  is 
completely  removed  as  shown  in  Figure  6d.  After  72  minutes  the 
initial  features  of  the  impacted  surface  can  no  Imger  be  identified. 

Another  sequence  of  micrographs  is  provided  in 
Figure  7  which  shows  a  triplet  of  interacting  ring  fractures  where 
the  conical  portion  of  one  ring  has  been  removed  by  crushing. 

This  region  was  first  observed  at  an  exposure  time  of  16  minutes. 

The  volume  of  material  removed  is  6  mils  by  6  mils  and  UO  microns  deep. 
There  was  no  change  in  the  configuration  at  2U  minutes  except  that 
the  fine- scale  chipping  has  been  active  and  is  beginning  to 
decorate  the  ring  fractures  already  present.  The  micrograph  of 
the  same  area  at  32  minutes  indicates  that  no  new  ring  fractures 
have  been  fonned  in  this  region.  The  micrograph  in  Figure  7  shows 
that  there  is  a  lot  of  loose  material  which  can  be  removed  easily 
on  subsequent  impacts.  It  is  at  this  stage  that  the  rate  of 
material  removal  should  be  uniform  and  proportional  to  the  number 
of  impacting  particles. 
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(a )  24  MINUTES 


(b)  32  MINUTES 


FIGURE  7  EROSION  IN  TOE  VICINITY  OF  A  LARGE  6  MIL 
PIT  ON  SPECIMEN  G1  (IlOX) 
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(c)  40  MINUTES 


(d)  56  MINUTES 

FIGURE  7  EROSION  IN  THE  VICINITY  OF  A  LARGE  6-MIL 
PIT  ON  SPECIMEN  G1  (IlOX) 
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An  example  of  a  spalled  region  Is  shown  in  Figure  8. 
Larger  areas,  approximately  10  mils  in  diameter,  were  occasionally 
removed  from  the  surface  in  the  early  stages  of  the  erosion 
process,  however  this  was  not  considered  a  major  contributor  to 
the  overall  weight  loss.  There  is  no  evidence  of  large-scale 
growth  at  this  nucleation  site,  but  a  fine-scale  chipping  away  of 
the  material  within  the  original  pit  and  in  the  surrounding  region. 
The  angularity  of  the  cleaved  surfaces  is  nearly  eliminated  after 
an  additional  exposure  of  24  minutes. 

The  observations  made  with  respect  to  Figures  6  to  8 
are  universal  and  are  applicable  to  all  of  the  large-scale  areas 
of  material  removal.  The  removal  of  the  interstitial  material 
between  adjacent  ring  fractures  is  negligible  and  is  no  longer  an 
operative  mechanism  once  the  relatively  short  exposure  time  is 
reeu:hcd  associated  with  a  stabilisation  of  the  ring  fracture  popu¬ 
lation.  No  large-scale  growth  was  observed  at  any  of  the  nucleation 
sites  where  measurable  material  removal  took  place  by  ring  fractures 
or  spallation.  The  original  conjecture  as  to  the  principal  erosion 
mechanism  in  fused  silica  using  2.8  mil  glass  beads  was  shown  to 
be  incorrect.  The  initiation  and  continuation  of  the  erosion 
process  is  due  entirely  to  the  fine-scale  chipping  associated  with 
the  2.8  mil  beads. 
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(a)  32  MINUTES 


The  fine-scale  chipping  process  can  be  seen  to  be 
preferential  and  until  the  time  the  original  surface  has  been 
removed,  somewhat  dependent  on  the  distribution  of  3  mil  ring 
fractures  on  the  surface.  Definite  growth  of  the  nucleated  areas 
can  be  seen  in  Figures  6  and  7.  This  process  is  in  accordance 
with  the  proposed  statistical  model.  However  there  are  certain 
limitations  w.ilch  make  it  extremely  difficult  to  evaluate  the 
data  required  in  the  proposed  theory  using  direct  microscopic 
observation  of  the  eroded  surface.  It  is  believed  that  these 
difficulties  can  be  overcome  in  the  future. 

Closer  examination  was  made  of  the  damage  associated 
with  what  has  been  referred  to  as  fine-scale  chipping.  At  higher 
magnification,  280X,  it  was  seen  that  the  small  pits  in  the 
range  of  slightly  less  than  0.1  rail  to  approximately  1  mil  were  in 
essence  material  removal  on  the  periphery  of  incomplete  ring 
fractures  whose  average  diameter  was  0.9  mils.  Complete  ring 
fractures  were  also  found,  but  usually  only  a  portion  of  the 
material  along  the  clrc\miference  is  removed.  These  fractures  were 
single  rings,  so  the  force  Imparted  by  the  particles  to  the  surface 
must  have  been  Just  sufficient  to  create  the  ring  fractures.  The 
general  character  of  the  damage  on  this  scale  is  shown  in  Figure  9. 
The  figure  shows  the  neighborhood  of  the  3-mll  rings  and  inter¬ 
actions  of  the  small  rings.  The  rings  are  formed  perferentlally 
and  seem  to  congregate  in  the  neighborhood  of  existing  rings.  The 


depth  of  the  materials  removed  is  only  3  to  10  microns  per  Impact. 
It  was  not  possible  to  count  all  of  the  pits  of  this  type  or  to 
record  their  locations  and  dimensions.  The  manner  of  growth  Is 
by  ring  fracture  Intersections,  but  alternative  methods  will  have 
to  be  devised  for  evaluating  the  area  removed  as  a  function  of 
time. 


The  weight-loss  data  for  specimen  G1  Is  given  In 
Figure  10.  After  an  exposure  time  of  70  minutes  the  rate  of  material 
removal  Is  reaching  a  steady-state  level.  It  can  also  be  seen 
that  only  3  mg  Is  removed  from  the  surface  at  the  exposure  time 
when  the  probability  for  the  occurrence  of  a  3  mil  ring  Is  zero 
(Figure  5).  This  provides  further  support  for  the  conclusion  that 
the  3  mil  ring  fractures  have  an  Insignificant  role  In  the  over- 
sdl  erosion  process. 

Specimen  02  was  exposed  to  the  11  mil  glass  beads 

i 

at  a  velocity  of  200  fps  smd  a  bead  flow  rate  of  0.016?  Ib/sec. 

The  saunple  will  experience  a  maximum  of  740partlcle  Impacts  per 
revolution  of  the  rotating  arm  according  to  Eq.  (22).  The  arm 
rotates  at  212  RPM  to  attain  a  linear  velocity  of  200  fps  at 
the  center  of  the  specimen.  Specimen  03  was  Impacted  by  the  11 
mil  beads  at  300  fps  (or  316  RFM)  at  the  same  bead  flow  rate  of 
0.0167  lb /sec.  The  erosion  process  will  be  described  for  both 
specimens  since  the  erosion  mechanisms  were  the  same  differing 
only  In  the  quantitative  aspects  of  the  process. 
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The  erosion  process  was  greatly  accelerated  when  the 
11  mil  beads  replaced  the  2.8  mil  beads  as  evidenced  in  the 
weight-loss  data  in  Figures  11  and  12.  Although  the  weight- los^t 
data  are  based  on  only  one  test,  the  characteristic  shape  can  be 
correlated  with  the  erosion  mechanisms  which  will  be  described  in 
a  chronological  manner.  The  erosion  mechanisms  operative  before 
reaching  the  uniform  erosion  rates  indicated  in  the  weight- loss 
data  are  quite  distinct  from  the  fine-scale  chipping  mechanism 
associated  with  the  2.8  mil  beads.  The  erosion  mechanisms  in  the 
early  stages  of  the  erosion  process  using  11  mil  beatds  are 
associated  primarily  with  conical  ring  fractures  and  subsequent 
impacts  on  the  damaged  surface. 

In  order  to  observe  a  relatively  low  density  of  ring 
fractures  on  the  surface  of  each  specimen,  they  were  exposed  to 
the  eroding  environment  for  a  time  increment  of  approximately  one 
second.  Surface  replicas  were  taken,  and  the  surfaces  were 
observed  optically  using  Noroarsky  phase  Interference.  The 
optical  micrographs,  Figures  13  and  14,  show  that  the  central 
region  of  the  inner  ring  is  not  damaged  by  the  impacts.  The 
circularity  of  the  ring  fractures  indicates  that  normal  Impacts 
are  taking  place.  It  was  not  possible  to  see  the  cone  fractures 
extending  into  the  interior  of  the  specimens  at  this  point; 
however,  subsurface  interactions  between  the  transient  stress 
fields  produced  by  particles  impacting  in  the  vicinity  of  already 
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Fijturc  1 2.  Wci|^|.Lm%  Dsla 
for  Specimen  G3 
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FIGURE  RING  FRACTURE  ON  THE  SURFACE  OF  SPECIMEN 

G3  AFTER  IMPACT  WIIH  11  MIL  BEADS  FOR  1  SEC. 

(125X) 
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existing  fracture  surfaces  could  be  observed  by  the  reflection  of 
light  from  these  surfaces.  A  good  example  of  this  effect  is  seen 
in  Figure  l4,  where  a  number  of  interactions  have  taken  place. 

The  ring  fractures  at  3OO  fps  have  outer  rings  which  are  approxi¬ 
mately  30  percent  larger  than  the  outer  rings  for  the  ring  fractures 
at  200  fps.  The  average  diameter  of  the  inner  rings  for  the 
fractures  formed  at  200  fps  and  3OO  fps  is  approximately  2.3  mils; 
however, the  range  of  the  inner  ring  sises  at  300  fps  is  fairly 
broad  (1.8  to  3.3  mils),  whereas  the  range  of  inner  ring  sires  at 
200  fps  is  relatively  narrow.  According  to  Hertr  theory,  Table  III, 
the  diameter  of  the  contact  area  should  be  2.86  mils  at  200  fps  and 
3.28  mils  at  300  fps.  The  diameter  of  the  outer  rings  in  each 
case  are  greater  than  the  calculated  values  for  the  contact  area. 

The  Rertrian  predictions  are  in  the  range  of  the  experimentally- 
determined  ring  fracture  diameters. 

Transmission  electron  micrographs  of  the  ring 
fractures  formed  in  the  early  stages  of  the  erosion  process  are 
shown  in  Figures  15  to  I8.  Since  the  stress  components  are  sym¬ 
metrical  about  the  contact  area,  the  initial  cracks  will  tend  to 
be  arcs  of  circles.  However,  the  relaxation  brought  about  by 
these  fractures  will  produce  a  redistribution  of  stress  with  the 
tendency  for  the  crack  to  propagate  in  a  direction  perpendicular 
to  local  trajectories  of  the  maximum  tensile  stress.  If  the 
loading  is  continued,  circular  fractures  will  be  produced  at 
other  flaws  In  the  same  manner  as  these  initial  cracks.  According 
to  Andrews^^^^,  if  the  crack  Is  able  to  follow  the  maximum  stress, 
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FIGURE  15  RING  FRACTURE  ON  SPECIMEN  G2  AFTER  AN  EXPOSURE  TIME  OF 
1  SECOND  (2000X) 
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it  will  not  be  circular.  If,  however,  for  any  reason  It  is 
prevented  from  doing  so,  it  may  follow  a  more  clearly  circular 
outline  for  a  short  time  and  then  branch  off  so  as  to  bridge  the 
distance  between  its  present  position  and  that  of  another 
circular  crack.  If  fracture  has  begun  at  several  places, 
development  of  each  crack  In  this  manner  is  to  be  expected 
giving  rise  to  the  complicated  figures  found  in  practice.  The 
micrographs  in  Figures  15  to  l8  are  indicative  of  this  process. 
Outside  of  this  general  observation,  no  particular  characteristic 
of  the  ring  fractures  investigated  using  the  electron  microscope 
has  been  identified. 

The  observations  of  the  specimen  surfaces  using 
Nomarsky  phase  interference  and  surface  replication  do  not  pro¬ 
vide  information  pertaining  to  the  characteristic  subsurface 
conical  fracture.  The  extent  of  the  conical  fracture  surfaces 
extending  into  the  interior  of  the  specimen  could  not  be  seen 
around  isolated  ring  fractures  occurring  in  the  very  early  stages 
of  the  erosion  process.  The  features  of  the  subsurface  fracture 
surface  did  become  evident  when  a  second  ring  fracture  occurred 
in  the  vicinity  of  an  existing  one.  Figure  19  shows  two  ring 
fractures  separated  by  a  distance  of  11-12  mils  with  no  obvious 
stress  wave  interaction  between  them.  Interactions  between  the 
dynamic  stress  field  of  an  Impacting  particle  which  results  in  a 
new  cone  fracture  and  the  flawed  surface  of  a  pre-existing  ring 
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(a)  NO  INTERACTION  BETWEEN  FRACTURE  SURFACES  i 


(b)  INTERACTING  SUBSURFACE  FRACTURES 

FIGURE  19  RING  FRACTURES  IN  PYREX  IMPACTED  BY  11  MIL  GLASS 
BEADS  AT  300  fps  (IlOX)  , 
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fracture  were  rarely  found  at  separation  distances  greater  than 
12  mils.  On  the  other  hand,  Figure  19b  shows  ring  fractures 
produced  under  the  same  conditions  separated  by  approximately 
8  mils  with  a  segment  of  the  common  fracture  surface  clearly 
visible.  This  effect  was  seen  repeatedly  In  the  microscopic 
observations.  In  nonoverlapping  rings  the  Interaction  generally 
produced  a  short  cracked  boundary  representing  the  Intersection 
of  the  conical  fracture  surface  associated  with  each  ring  fracture. 
This  clearly  visible  demarcation  Is  perpendicular  to  the  line 
between  the  centers  of  the  two  adjacent  ring  fractures.  A 
particle  Impacting  Just  outside  the  periphery  of  a  ring  fracture 
produces  a  slightly  different  effect  as  shown  In  Figure  20.  Over¬ 
lapping  ring  fractures  show  surprisingly  little  Influence  on 
each  other  as  seen  In  the  Nomarsky  interference  micrographs, 
Figures  13  and  l4,  and  in  Figure  20.  As  more  particles  impact 
the  neighborhood  of  the  pre-existing  ring  fractures  a  complex 
pattern  of  subsurface  fractures  develop  as  shown  in  Figure  21. 

The  micrograph  In  Figure  22a  is  taken  at  the  surface  of  the  spec¬ 
imen  while  the  micrograph  in  Figure  22b  is  focused  at  a  point 
below  the  surface  which  shows  the  straight-line  segments  produced 
In  the  manner  described  above.  The  suDsurface  fracture  process 
continues  without  any  significant  material  removal  from  the 
surface.  Subsequent  Impacts  on  the  interstitial  regions  between 
cone  fractures  produce  an  increasingly  complex  network  of  fracture 
surfaces,  but  no  new  cone  fractures  are  produced  in  the  highly- 
flawed  interstitial  areas.  The  general  appearance  of  the  surface 


-74- 


PTaimE  20  OVERLAPPLING  AND  NEARLY  OVERLAPPING  RINGS 
FRACTURE  IN  PYREX  IMPACTED  BY  11  MIL  BEADS 
AT  400  fps  (IlOX) 


FIGURE  21 


SUBSURFACE  FRACTURE  INTERACTIONS  IN  PYREX 
IMPACTED  BY  11  MIL  BEADS  AT  400  fps  (90X) 
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(a)  FOCUSED  AT  SURFACE 


(b)  FOCUSED  BELOW  SURFACE  SHOWING  FRACTURE  SURFACE 
INTERSECTIONS 

FIGURE  22  SUBSURFACE  FRACTURE  INTERACTION  IN  PYREX  IMPACTED 
BY  II  MIL  BEADS  AT  400  fps  (90X) 
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of  specimen  G2  at  15  seconds  and  60  seconds  Is  Indicated  in 
Figure  23.  The  small  white  spots  on  the  micrograph  are  evidence 
of  collisions  which  did  not  result  In  ring  fractures.  The 
particle  bombardment  has  relatively  little  Influence  on  the 
Interior  surface  region  of  the  ring  fractures.  The  interstitial 
regions  after  an  exposure  time  of  60  seconds  are  still  intact 
over  most  of  the  surface  area  of  the  specimens. 

The  distribution  of  diameters  of  the  inner  rings  of 
the  ring  fracture  produced  in  specimens  G2  and  G3  is  shown  in 
Figures  24  and  25,  respectively.  These  data  were  obtained  at  an 
intermediate  stage  of  the  overall  exposure  for  each  specimen  and 
are  based  on  measurements  of  200  rings  on  the  surface  of  each 
specimen.  The  average 'diameter  for  the  inner  rings  for  specimen 
G2  is  2.63  mils,  while  it  is  2.51  mils  for  specimen  G3.  This 
result  is  consistent  with  the  average  diameter  of  the  inner  rings 
evaluated  in  the  very  early  stage  of  the  erosion  process.  The 
reason  for  the  average  ring  size  being  larger  at  the  lower  Impact 
velocity  is  believed  to  be  due  to  measuring  errors. 

Surveying  the  central  region  of  specimen  G2,  it  was 
found  that  the  density  of  ring  fractures  per  square  inch  was 
3336  at  15  seconds,  4l40  at  60  seconds  and  2930  at  I30  seconds. 
Similar  data  for  specimen  G3  show  that  the  ring  fracture  density 
(rings /square  inch)  is  I820  at  3  seconds  and  i960  at  18  seconds 
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reaching  a  maximum  value  within  the  given  time  range.  The 
probability  for  the  formation  of  ring  fractures  is  quite  high 
under  these  impact  conditions.  From  these  data  it  is  seen  that 
the  maximum  number  of  ring  fractures  is  established  quite  early 
in  the  erosion  process  and  that  a  majority  of  the  rings  formed 
survive  the  subsequent  bombardment  of  particles  as  the  erosion 
process  continues.  The  reasons  for  the  difference  in  the  density 
of  cone  fractures  between  the  two  specimens  and  the  removal  of 
certain  cones  have  been  established  from  the  experimental 
observations.  The  explanation  is  illustrated  in  Figure  26. 

The  top  surface  appears  as  in  Figure  26a.  However,  a  cross 
section  of  the  conical  fracture  surfaces  indicates  that  cones  A 
and  C  intersect  below  cone  B.  Such  conditions  have  been  observed 
microscopically.  Cone  B  is  only  partially  attached  to  the  bulk 
material  while  cones  A  and  C  are  firmly  anchored.  The  colliding 
particles  are  capable  of  producing  fractures  around  cone  B 
while  weakening  its  existing  support  and  it  ultimately  breaks 
loose.  The  depth  D  is  found  to  be  fairly  constant  over  the 
entire  eroded  area.  In  specimen  G2  D  is  approximately  50-60  microns 
while  in  specimen  G3  the  depth  D  is  generally  80  to  90  microns. 

The  larger  D  in  specimen  G3  is  attributed  to  the  higher  velocity 
of  approach.  It  is  the  greater  depth  of  fracture  penetration 
for  specimen  G3  which  contributes  to  the  marked  difference  in 
ring  fracture  density  when  compared  with  specimen  G2. 
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(a)  RING  FRACTURES  ON  GLASS  SURFACE 


FIGURE  26  MECHANISM  FOR  REMOVAL  OF  CONE  FRACTURES 
FROM  SURFACE 


Following  the  erosion  process  in  specimen  G2  it  is 
seen  that  around  an  exposure  time  of  90  seconds  large-scale 
removal  of  material  from  the  interstitial  regions  between  the 
firmly- anchored  truncated  cones  is  the  primary  erosion  mechanism. 
The  glass  in  the  interstitial  areas  is  severely  fractured  and 
sizable  areas  are  removed  during  each  time  Increment.  The  micro¬ 
graph  in  Figure  27  Illustrates  this  process.  Just  above  the 
center  of  Figure  27  there  is  a  triplet  of  cones  from  which  the 
interstitial  glass  has  been  removed.  In  the  vicinity  of  the 
triplet  the  Imprint  of  two  cones  which  did  not  survive  the 
impact  loadings  are  also  clearly  visible.  The  conditions  for 
their  disintegration  correspond  to  that  shown  schematically  in 
Figure  26.  A  number  of  other  cones  which  are  going  td*be 
removed  with  the  interstitial  glass  can  alsij  be  seen  in  Figure  27. 
The  sequence  of  micrographs  in  Figure  28  shows  the  growth  of  an 
eroded  area.  The  cones  protruding  from  the  bulk  material  are  not 
always  visible  in  the  micrographs;  only  the  interstitial  material 
is  being  removed,  the  cones  at  this  stage  are  relatively  unaffected 
by  the  particle  impacts.  The  glass  comprising  the  truncated  cones 
is  quite  strong  and  fairly  resistant  to  erosion.  A  number  of 
areas  on  the  sample  surface  were  observed  microscopically  during 
this  stage  of  the  erosion  process,  and  it  was  found  that  the 
growth  behavior  of  nucleated  regions  was  quite  irregular.  At 
130  seconds  most  of  the  interstitial  material  was  removed  from 
the  primary  erosion  area  on  the  specimen  surface  as  shown  in 
Figure  29.  The  removal  of  the  interstitial  material  completes 
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FIGURE  27  REMOVAL  OF  FRACTURED  GLASS  FROM  INTERSTITIAL 
REGIONS  ON  SURFACE  OF  SPECIMEN  G2  (IlOX) 


(a)  90  SECONDS 
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FIGURE  28  GROWTH  OF  AREA  OF  INTERSTITIAL  GLASS  REMOVAL 
ON  SURFACE  OF  SPECIMEN  G2  (IlOX) 
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(c  )  110  SECONDS 
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FIGURE  28  GROWTH  OF  AREA  OF  INTERSTITIAL  GLASS  REMOVAL  ON 
SURFACE  OF  SPECIMEN  G2  (IlOX) 
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another  phase  of  the  erosion  trocess  which  overlaps  with  the  final 
stage. 


The  general  appearance  of  the  eroded  surface  Is  now 
comparable  to  a  waffle  Iron.  The  appearance  of  the  surface  of 
each  specimen  Is  shown  at  lower  magnification  In  Figure  30.  Solid 
cones  are  protruding  from  the  surface.  The  Impacting  particles 
do  little  damage  In  a  single  Impact.  The  sides  of  the  cones  and 
common  points  on  Interaction  are  highly  resistant  to  material 
removal.  The  particles  must  pulverize  the  glass  In  the  cones 
beginning  at  the  plateau  and  slowly  work  down  to  the  bulk  material 
(Figure  31)*  As  seen  In  the  plots  of  the  weight-loss  data 
(Figures  11  and  12)  this  process  continues  at  a  uniform  rate. 

No  new  erosion  mechanisms  come  Into  play. 

The  conditions  of  the  glass  beads  were  checked 
microscopically  after  colliding  with  the  specimens  In  the 
erosion  tests,  and  It  was  found  that  no  fragmentation  of  the 
beads  was  taking  place  and  that  very  few  fractures  could  be 
found  on  the  surface  of  the  beads. 
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(a)  SPECIMEN  G2  AT  100  SECONDS  (25X) 


(b)  SPECIMEN  G3  AT  6  SECONDS  (20X) 


FIGURE  30  GENERAL  APPEARANCE  OF  ERODED  SURFACE  OF 
SPECIMEN  G2  AND  G3 
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VI.  GENERAL  DISCUSSION 


In  order  to  couple  the  general  erosion  model  developed 
In  Section  III  with  the  erosion  behavior  of  specific  materials, 
an  experimental  procedure  was  devised  for  evaluating  the 
statistical  parameters  appearing  In  the  theory  for  the  special 
case  of  spherical  glass  beads  colliding  with  the  plane  surface 
of  a  specimen  of  fused  silica.  The  experimental  results  are  now 
available  for  the  future  development  of  the  erosion  model. 

It  was  pointed  out  In  Section  V-A  that  the  applicability 
of  tne  Hertzian  analysis  may  be  questioned  on  the  basis  of  the 
assumptions  that  the  contact  area  be  small  In  comparison  to  the 
diameter  of  the  Impinging  sphere  and  that  the  duration  of  Impact 
be  long  In  comparison  to  the  time  It  takes  stress  waves  to  traverse 
the  Impacting  bodies.  However  the  Hertzian  prediction  of  the 
size  of  the  contact  area  Is  found  to  be  In  reasonable  agreement 
with  the  experimental  observations  pertaining  to  the  production 
of  ring  fractures  on  the  glass  surfaces.  In  specimen  G1  a  high 
probability  exists  for  the  formation  of  single  ring  fractures 
approximately  0.9  mils  in  diameter.  According  to  the  Hertz  theory 
the  maximum  diameter  of  the  contact  area  will  be  O.70  mils. 

In  static  loading  of  Indenters  on  glass  plates  it  Is  found  that 
ring  fractures  generally  occur  outside  of  the  contact  area  by 
diameters  up  to  50  percent^  greater  than  the  contact  diameter. 
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The  usual  difference  is  10  to  20  percent.  The  experimental  work 
of  Tsai  and  Kolsky^^^^  on  Impact  of  steel  ball  bearings  with 
radii  ranging  from  1  '32  to  1/4  inch  on  plate  glass  Indicates 
that  the  fracture  stress  in  the  dynamic  tests  is  approximately 
one- third  higher  than  the  static  fracture  stress.  These  obser¬ 
vations  tend  to  make  the  calculated  radius  of  the  contact  area 
more  closely  approximate  the  experimentally-evaluated  radius 
for  the  ring  fractures. 

The  Hertz  theory  predicts  the  diameter  of  the  contact 
area  for  specimens  G2  and  03  should  be  about  5  mils.  A  very 
high  probability  was  found  for  the  occurrence  of  multiple  ring 
fractures.  The  average  diameter  of  the  innermost  ring  was 
approximately  2.5  mils:  the  diameter  of  the  outer  rings  in  each 
case  was  determined  fran  the  experiments  to  be  slightly  larger 
than  the  calculated  value  for  the  diameter  of  the  contact  area. 
The  comparison  between  theory  and  experiment  is  reasonable. 

The  stress  generated  at  an  intermediate  stage  of  the  compression 
phase  of  the  impact  is  sufficient  to  initiate  ring  fractures,  but 
the  magnitude  of  the  load  continues  to  act  along  with  the  expan¬ 
sion  of  the  contact  area.  Additional  ring  fractures  are  produced 
as  a  larger  area  reaches  the  critical  stress  level.  The  process 
stops  when  the  maximum  contact  area  is  reached  where  it  is  still 
possible  for  a  crack  to  develop  Just  outside  this  region  before 
the  magnitude  of  the  load  begins  to  subside.  The  static  experi¬ 
ment  of  Sheldon  and  Finnie^®^  for  the  variation  of  the  center  of 


approach  of  a  1/8  inch  steel  ball  on  a  glass  plate  indicates 
that  Hertzian  expression  for  the  approach  distance,  Eq,  (l4) 
deviates  from  the  experimental  value  by  only  25  percent  when  the 
applied  load  is  more  than  120  times  the  load  for  the  initial 
ring  fracture.  It  thus  appears  that  cracked  material  will  trans¬ 
mit  load.  Sheldon  and  Flnnle  state  that  this  is  due  to  the  fact 
that  the  stresses  in  the  contact  region  are  primarily  compres¬ 
sive. 


Only  Indirect  evidence  can  be  provided  in  regard  to  the 
second  assumption  of  the  Hertz  theory.  According  to  the  calcu¬ 
lation  of  the  dynamic  contact  stresses  at  the  boundary  of  the 
area  of  contact  provided  by  Tsal^^^\  it  is  seen  that  the  radial 
tensile  stress  component  is  essentially  equal  to  the  Hertzian 
radial  tensile  stress  for  a  maximum  radius  of  contact  less  that 
10  mils  Independent  of  the  duration  of  contact.  The  size  of  the 
contact  areas  found  experimentally  are  considerably  smaller  than 
this  value. 

The  measurements  of  the  ring  fractures  would  tend  to 
support  the  Hertzian  evaluation  of  the  contact  stresses  for  the 
Impact  conditions  appropriate  to  the  erosion  tests  conducted 
in  this  report.  The  Hertzian  expression  for  the  contact  stresses 
will  therefore  be  used  to  construct  an  appropriate  fracture 
criterion  for  the  small  bead  sizes  used  in  the  experimental 
program.  This  task  will  be  considered  in  future  work. 
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Two  additional  observations  can  be  made  pertaining  to 
ring  fractures  based  on  these  initial  studies.  First,  the 
semi-angle  of  the  conical  fracture  surfaces  exposed  in  the 
intermediate  stage  of  the  erosion  of  specimens  G2  and  G3 
appeared  to  be  closer  to  45°  rather  than  65.5°  in  fused  silica 
as  measured  by  Benbow^^^^.  Secondly,  the  anomolous  occurrence 
of  3  nill  ring  fractures  on  the  surface  of  specimen  G1  cannot 
be  explained.  Their  existence  is  not  at  all  in  conformity  with 
the  Hertzian  theory.  The  3  mil  ring  fractures  were  incorrectly 
taken  to  be  the  principal  element  contributing  to  the  nucleatlon 
of  large-scale  pit  growth  on  the  surface.  The  number  and  config¬ 
uration  of  the  3  niil  ring  fractures  was  tabluated  after  each 
erosion  Increment,  until  it  was  discovered  that  an  exposure 
time  was  reached  when  no  new  3  mil  ring  fractures  were  being 
formed.  It  then  became  evident  that  the  smaller  ring  fractures 
described  above  were  the  primary  mechanism  of  material  removal. 

Very  small  pieces  of  material  were  removed  along  the 
periphery  of  the  0.9  mil  rings.  Once  a  region  began  to  develop, 
these  rings  overlapped  each  other  and  exhibited  the  form  of 
preferential  growth  described  by  the  theory.  They  were  found  to 
cluster  in  the  neighborhood  of  the  3  mil  rings.  The  complex 
spatial  growth  pattern  was  difficult  to  tabulate  by  the 
procedures  established  for  the  anticipated  erosion  damage  on  a 
somewhat  larger  scale.  The  depth  of  the  surface  layer  removed 
was  only  five  to  ten  microns,  so  the  erosion  of  layers  below  the 


surface  layer  was  nearly  impossible  to  follow.  The  chipping 
associated  with  the  ring  fracture  was  responsible  for  material 
removal  for  almost  the  full  duration  of  the  cummulative  erosion 
time.  Ring  fractures  were  still  being  formed  on  the  eroded 
surface  Just  prior  to  the  attainment  of  the  steady- state  erosion 
rate  (Figure  10).  When  the  steady-state  erosion  rate  was  attained, 
material  removal  took  place  by  fine-scale  chipping  of  pieces 
from  the  highly-flawed  surface  without  any  evidence  of  the  ring 
formations . 

Flnnle^^^^  developed  a  model  for  the  erosion  of  brittle 

materials  based  on  the  idea  that  the  volume  removed  per  spherical 

bead  Impact  V  ,  Is  proportional  to  the  product  of  the  statistical 
^  - 

average  of  the  square  of  the  outermost  ring  crack  (a*^),  and 
the  depth  of  the  Inner  ring  fracture.  Using  the  statistical 
flaw  theory  of  Oh  and  Flnnle'  ’  and  the  usual  Hertzian  relations, 
Plnnle  finds 

Vp  OC  „  P  (24) 

The  experimental  confirmation  of  this  result  based  on  the  erosion 
of  glass  Indicates  that  the  linear  proportionality  of  Vp  with  F, 
the  grouping  of  parameters  In  Eq.  {2h),  was  not  upheld,  although 
the  data  did  cluster  along  a  straight  line  with  a  slope  of  I.30 
on  a  log- log  plot  of  Vp  vs  F. 
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If  Vp  is  evaluated  by  the  conventl<nal  solid  particle 
erosion  tests,  it  would  presumaoly  correspond  to  the  steady-state 
rate  of  erosion.  Eq.  (24)  in  essence  relates  conditions  in  tl  e 
very  early  stages  of  the  erosion  process  to  measurements  of  the 
volume  removed  (or  equivalently,  weight  loss)  in  the  later  sta 
of  the  erosion  process.  According  to  the  erosion  mechanism 
studies  for  fused  silica,  it  is  seen  that  Finnic' s  model  is 
representative  of  the  primary  mechanism  for  material  removal  in 
specimen  G1  over  the  entire  period  of  exposure.  However,  the 
erosion  of  specimen  02  and  03  Involves  several  erosion  mechanisms. 

The  prevalent  mechanism,  chipping  of  the  conical  frustum  popula¬ 
tion,  for  the  steady-state  region  of  the  weight  loss  data  (Flnires 
11  and  12)  is  entirely  different  from  the  ring  fracture  mechanism 
which  is  the  basis  of  Finnie's  model.  In  this  case  agreement 
between  the  experimental  and  calculated  values  of  would  be 
fortuitous. 

These  initial  experiments  in  this  report  have  provided  valuable 
insights  into  the  role  of  particle  size  and  velocity  on  the 
erosion  characteristics  of  Inorganic  glasses.  These  experimental 
observations  will  be  used  in  conjunctl  n  with  the  available 
Hertzian  analyses  of  Impact  phenomena  to  further  advance  the 
quantitative  theory  of  erosion  of  brittle  materials. 
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VII.  CONCLUSIONS 


The  initial  studies  on  analytical  modeling  of  solid 
particle  erosion  have  demonstrated  that  it  is  feasible  to  Include 
more  details  of  the  erosion  behavior  of  materials  in  the  general 
statistical  framework  which  is  developed  in  this  report.  This 
Is  possible  through  simplification  of  the  Impact  conditions, 
direct  microscopic  observation  of  the  time  evolution  of  erosion 
damage  on  material  surfaces,  and  the  preliminary  observation 
that  the  Hertzian  theory  may  be  applicable  to  impact  conditions 
on  a  scale  appropriate  to  solid  particle  erosion. 

At  the  present  time  there  Is  a  general  lack  of  experimental 
results  pertaining  to  the  critical  velocity  for  fracture  of 
Inorganic  glasses  and  other  brittle  materials  when  Impacted  by 
spherical  particles  less  than  2k  mils  in  diameter.  This  size  of 
Impacting  particle  Is  approaching  the  upper  limit  of  particle 
sizes  which  are  of  Interest  In  Investigation  of  sand  erosion. 

Tests  are  being  planned  for  evaluating  the  most  probable  value  of 
the  critical  velocity  for  fracture  as  a  function  of  particle 
sizes  In  the  range  from  1  to  l8  mils  using  the  APML-Bell  erosion 
apparatus  for  the  materials  which  are  selected  for  study  In  this 
research.  The  quantitative  evaluation  of  the  critical  velocity 
Is  required  for  obtaining  analytical  representations  of  the 
primary  erosion  mechanisms  which  can  be  Incorporated  Into  the 
functional  forms  of  the  pit  nucleatlon  and  growth  functions. 
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In  conjunction  with  the  detailed  studies  of  direct  Impacts, 
a  series  of  controlled  experiments  are  planned  to  progressively 
increase  the  complexity  of  the  loading  conditions  using  the  AFML- 
Bell  erosion  apparatus.  This  effort  is  Intended  to  extend  the 
range  of  application  of  the  theoretical  predictions  to  conditions 
which  are  of  greater  Interest  in  solid  particle  erosion.  The 
impact  conditions  to  be  considered  are  oblique  impacts,  impacts 
by  arbitrarily- shaped  particles,  higher  velocity  impacts,  and 
finally  erosion  by  liquid  drops.  More  complex  forms  of  material 
behavior  than  that  of  the  inorganic  glasses  will  also  be  investi¬ 
gated  once  viable  predictive  tools  are  available. 
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APPENDIX 


COMPUTER  PROGRAM  FOR  ANALYTIC  MODELING 

The  system  of  equations  governing  the  erosion  of  solid 
surfaces  by  the  nucleatlon  and  growth  of  erosion  pits,  Eq.  (11), 
are  programmed  for  digital  computer  computations.  The  specific 
input  data  for  the  program  provided  is  highly  restricted  in 
order  to  make  It  explicit.  There  Is  no  difficulty  In  including 
a  very  broad  range  of  Input  functions  now  that  means  are  available 
for  evaluating  the  Integrals  In  Eq.  (11).  The  form  of  the  output 
from  the  computer  computations  is  also  arbitrary. 

Generalized  forms  of  the  program  would  Include  a  different 
nucleatlon  probability  function  and  growth  rate  function  for  each 
of  the  eroding  layers  whose  thickness  can  be  varied  as  material  Is 
removed  from  the  bulk  solid.  Since  It  Is  our  purpose  to  specify 
this  information  based  on  experimental  and  theoretical  considera¬ 
tions  appropriate  to  realistic  material  behavior  and  Impact 
conditions,  relatively  simple  forms  of  the  input  data  were  selected 
to  check  If  the  program  was  operational.  This  has  been  verified 
assuming  the  relevant  functions  are  identical  for  all  layers  of  the 
laminate  model  which  are  all  the  same  thickness.  Identical  linear 
growth  rates  are  assumed  In  both  lateral  directions  on  the  surface 
of  the  eroding  solid,  and  the  pit  nucleatlon  function  Is  simply 
proportional  to  the  exposure  time. 
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ALX?  xSL'MUHl*  IY(I-2)  £AL>1  (Y(I  )) 
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4  I6R«-l 
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